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Abstract 

The pomeron-pomeron interaction is studied in the perturbative approach of BFKL- 
Bartels. The total pomeron-pomeron cross-section is proportional to a* / yjt\t2 where ^fs 
is the cm. energy and t± t 2 are the virtualities of the colliding pomerons. Upon calculating 
the coefficient the cross-section is found to be of the order 2.2 mb at \/s = 6 TeV and 
y/—t{ ~ ~ 1 GeV/c. 
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1 Introduction 

The study of high-mass diffractive events brings into consideration properties of the pomeron 
(P), which may be pragmatically associated with a high-mass object appearing between the 
rapidity gaps. In particular, diffractive events in DIS allow to study the structure function of 
the pomeron and central diffractive events in hadronic or photonic collisions reveal the prop- 
erties of the PP interaction. Most of the theoretical activity in this domain has concentrated 
on the problem of extracting the pomeron properties from the data, rather than to study 
these properties themselves. 

Such a study is possible in the perturbative framework of the "hard pomeron", in which 
the realistic QCD is modified by introducing an infrared cutoff (e.g. via the Higgs mechanism) 
and subsequently fixing a (small) coupling constant. As is well known, this theory results 
infrared stable in the colourless sector, so that one might think that at least part of its 
predictions remain valid for the realistic QCD. 

In a recent publication [1], in this approach, we have studied the triple pomeron interac- 
tion, which is vital for the pomeron structure function. The present note extends this study 
to a more complicated case of the PP scattering cross-section, which evidently involves two 
triple pomeron vertices. 

Since the hard pomeron theory is scale invariant, the dependence of the PP scattering 
cross-section on the relevant variables is trivially predicted in this theory. If the two colliding 
pomerons have their momenta I and I' with t\ = I 2 < and t<i = I' 2 < 0, and their cm. 
energy squared is s = (I + l') 2 , then from dimensional considerations and properties of the 
BFKL pomeron [2] one finds for the total PP cross-section 

« 4 s A 

a%(s, tl ,t 2 ) = c s . (1) 

\Jat\ti Ins 

Here the pomeron intercept A and the parameter a are known functions of the strong coupling 
constant a s (see Eq. (18)) and c is an unknown number. Therefore the problem of studying 
the PP cross-section reduces to finding the numerical constant c in (1). One should have in 
mind that this constant factor is in fact rather ill-defined. On the one hand, the separation of 
the PP cross-section from the rest of the amplitude inevitably involves a certain arbitrariness, 
known as a choice of the "flux factor" [3]. In the hard pomeron theory, as we shall see, this 
problem is aggravated because the pomeron coupling to external sources results singular for 
small momentum transfers. On the other hand, the scale of In s in the hard pomeron approach 
remains arbitrary, which introduces an undetermined factor in (1). 

In the following we calculate the coefficient c within a certain (rather obvious) choice of 
the flux factor and neglecting the numerical factor which comes from s A in (1). This latter 



3 



approximation corresponds to the (theoretically justified) assumption of a very small A. In 
applications we take the scale involved of the order of 1 GeV 2 , as favoured by the pomeron 
phenomenology. Our calculations give a rather large value for c: 

c ~ 35370. (2) 

Possible experimental consequences of this result are discussed in Section 4 of this note. 
Sections 2 and 3 are dedicated to its derivation. 



2 The central diffraction cross-section 

The PP cross-section enters as a part of a more general central diffraction cross-section 
in which the colliding hadrons (or photons) produce a high mass M in the central region 
separated from the target and projectile by large rapidity gaps (Fig. 1). As seen from this 
figure, both projectile and target emit pomerons with momenta I and I', respectively, which 
interact in the central part of the diagram. We introduce the standard energetic variables 
(see Fig. 1 for the notation of the momenta involved) 

S=(P1+P2) 2 , S 12 = (p 1 + l') 2 , S 21 = (p 2 + l) 2 , M 2 = (/ + /') 2 , (3) 

constrained by the relation 

S12S21 = sM 2 . (4) 

We assume that all these variables are large with s » si 2 , S21 >> M 2 . We also introduce 

Sl = s 2 i/M 2 , s 2 = s 12 /M 2 , (5) 

which serve as energetic variables for the upper and lower pomerons in Fig. 1, respectively. 

In the hard pomeron approach the pomerons which connect the central part of Fig. 1 
with the target and projectile are identified as BFKL pomerons. As to the central part itself, 
it is also given by a BFKL pomeron coupled to the upper and lower pomerons by two Bartels 
vertices K 2 ^4, which describe transition of two reggeized gluons into four [4] (see Fig. 2). 
Explicitly the upper vertex K 2 ^4, apart from a colour factor and a factor g A , is given by 

#2-4 -qi; q 2 , 1 - <?2, -/ + q 3 , -93) = K i (q 1 ,q 2 , q 3 ) = 

gjg| + gigf qKq2 - q 3 ) 2 ^ 

(Q1-Q2) 2 {qi-qz) 2 {qi- q2) 2 {qi- q?,) 2 ' 

The lower vertex has the same form with the primed momenta. 

Then we obtain the expression for the absorptive part D corresponding to Fig. 1 in the 
form 

D = (l/16) 5 16 iV 4 (iV 2 - l)( S 2 /M 2 )(2vr)- 12 f f[(d 2 q t d 2 q' l )K l (q 1 ,q 2 ,q 3 )K l ,(q' 1 ,q 2 ,q> 3 ) 

J i=i 
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M*i, l, 92)^1 (si, -I, -q 3 )<h(s2,l', q' 2 )Ms2, -l', -q' 3 )G(M 2 ,0, Ql , -q[). (7) 

Here 4>K2) (s, I, q) is a pomeron (a solution to the BFKL equation) coupled to the projectile 
(target) colourless external source with energy squared s, total momentum I and one of 
the gluon's momentum q. The function G(M 2 , 0, q\, q[) is the BFKL Green function for the 
energy squared M 2 , total momentum zero and initial and final momentum of one of the gluons 
qi and q[, respectively. It is assumed that the colour factor for each source is (l/2)5 a & and 
that each source is proportional to g 2 . The factor (l/16)iV 4 (iV 2 — 1) comes from the colour 
variables, TV being the number of colours (N = 3 for the physical QCD). Our normalization 
is that the double inclusive cross-section described by Fig. 1 is given by 

d(T = D (8) 

dt 1 dt 2 ds 12 ds 2 i 256vr 4 s 3 ' 1 ; 

Since for I ^ the solution of the BFKL equation is easier to obtain in the (transversal) 
coordinate space, we pass to this space by presenting 

(f>(s, l,q)= J d 2 r<p(s, I, r) exp[ir(q - 1/2)]. (9) 

Evidently r is the transversal distance between the gluons. Then (7) transforms into 

D = (l/16)g ie N\N 2 - l)(s 2 /M 2 ) J {[(d^d^K^^r^K^r^r'^r',) 

i=l 

ex P HZ(r 2 + r 3 )/2 - il'(r' 2 + r' 3 )/2] 
<h(*i,h r 2 )Msi, -I, ~r 3 )Ms2, 1', r' 2 )Ms2, -I', -r' 3 )G(M 2 , 0, n, -r[), (10) 
where the 2 to 4 reggeon vertex in the coordinate space is obtained from (6) to be 

Ki(n,r 2 ,r 3 ) = -2 



5 2 (r 2 )5 2 (r 3 )V 2 5 2 (r) + ^ 2 (r3K 2 (r 2 Vi)V^ 2 (r) + 



^S 2 (r 2 )r^ 2 (r 3 V!)V 2 5 2 (r) + ^ r^ 2 (r 2 r 3 )V\5 2 (r) 



(11) 

with r = r\ + r 2 + r 3 . The solutions (f> vanish when r = 0. So only the last term in (11) 
survives. We then can rewrite (10) as 

D = (l/4)g w N 4 (N 2 - l)(s 2 /M 2 )(2ir)' 8 

J d 2 qd 2 q'x(M 2 , q + 1/2, q' + l> /2)[V qX i(si, I, q)] 2 [V ' q 'X2(s 2 , l', q')} 2 , (12) 

where 

X(s, q, q') = J d 2 rd 2 r exp(iqr + igV)V 4 V /4 G(M 2 , 0, r, -/) (13) 

and 

Xi, 2 {s,l,q)= d 2 rr~ 2 4>i )2 (s,l,r)exp(iqr). (14) 
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To calculate the function x( s iQil') we first integrate by parts in (13) to remove the 
derivatives, 

x (s, q, q') = {qq' f J d 2 rd 2 r'G{M 2 , 0, r, -r') exp(iqr + iq'r'). (15) 
The BFKL Green function at / = and large s is given by the expression [5] 

c (s ,o,,- r <) = ^ ^-^r- (16) 

where 

u{y) = (g 2 N/2TT 2 )[ip{l) - Re^(l/2 + iu)]. (17) 
Small values of v play the dominant role in (16) at large s, so that we can approximate 

u(y) = A- av 2 ; A = (g 2 N/i: 2 ) In 2, a = (7g 2 N/2tt 2 )((3). (18) 

Calculating the integrals over r and r' in (15) and then the remaining integral over v we 
obtain 

x(5 ,,,, = 2^ /^ exp (_!5!w)V (19) 

V a Ins V ams y 
At large s the exponential factor in (19) can evidently be neglected. 

Now we turn to the functions Xi,2(s,l,q)- The solutions 1( - 2 ) can be obtained by using 
the Green function of the BFKL equation for a given total momentum G(s, I, r, r'). Following 
[6] and taking into account that the Green function G vanishes if r or r' are equal to zero 
and is isotropic at high energies one finds for, say, the projectile: 

01 (s, I, r) = -2 J d 2 rG(s, I, r, r)pi{l, f) exp(ilr/2). (20) 

Here pi(l,r) is the colour density of the projectile as a function of the intergluon distance r 
with the colour factor (l/2)S a b and g 2 separated. The explicit form of p can easily be found 
if the projectile is a highly virtual photon. For the longitudinal photon then 

, n 4.^21 21 N f ,1 

Pi\r) = E Z f I ^H 1 " a)] 2 K 2 (e f r)eM-^lr), (21) 

where e 2 = \p 2 \a(l — a) + m 2 and rrif and Zf are the mass and charge of the quark of flavour 
/, respectively. For the transverse photon a slightly more complicated formula emerges 



2 N f rl 



[m)Kl{e f r) + [a 2 + (1 - a) 2 )e) K?(e/r) - a(l - a)(l - 2a)K (e / r)zZV r K (e / r)} exp(— iaZr) 

(22) 



6 



Eqs (21), (22) generalize the well-known formulas for zero momentum transfer [6]. For a 
hadron projectile or target the form of the density is, of course, unknown, although one 
expects that it should be less singular at r = and normalizable. 

The leading contribution to the BFKL Green function at I ^ has the form [5] 

1 f duu 2 

G^r') = W¥ j {v2 + 1/4)2 s^Er(r)Er(r% (23) 

where 

// r \ l+2iv 

rf2i?eXP( " R) U + r/2||fl-r/2| ) ' <24) 
It is evident that the Green function (23), transformed into momentum space, contains terms 
proportional to 5 2 (l/2±q), which should be absent in the physical solution (this circumstance 
was first noted by A.H.Mueller and W.-K.Tang [7]). For that, (23) goes to zero at r = 0. 
Terms proportional to 5 2 {l/2 + q) are not dangerous to us: they are killed by the vertex ^2^4 • 
To remove the dangerous singularity at q = 1/2 and simultaneously preserve good behaviour 
at r = we therefore make a subtraction in E, changing it to 

r T / r \ 1 + 2iv 

Et(r) = / [( |fl + r/2||fl _ r/2| ) - I* + r/2|-'— + I* - r/2|— - 

(25) 

This subtraction removes the 5 singularity at q = 1/2 and doubles the one at q = —1/2, the 
latter eliminated by the kernel i^2-^4- 

Integration over r in (14) leads to the integral 

J(l, q)= J (d 2 r/(27r) 2 )(l/r 2 )E?(r) exp(igr). (26) 

This integral is convergent at any values of v, the point v = included, when the convergence 
at large values of r and R is provided by the exponential factors. So in the limit s — > 00, 
when small values of v dominate, we can put v = in J: 

,(,„, = /(WW') r ,^^| + (1 /0.n («) . (27, 

The second term comes from the subtraction terms in (25). Passing to the Fourier transform 
of the function \jr one can represent (27) as an integral in the momentum space 

J(/ ' Q) ~ m J (2,)\l/2 +P \\l/2- P \\q +P \ ■ (28) 

In the same manner we can put v = in the function Ef u (f) obtaining for the integral 
over the transverse dimensions of the projectile 

f d 2 Rd 2 rexp[il(R + r/2)]r Pl (l,r) 



I 



\R + r/2\\R-r/2\ 



= irR 1 F 1 (t 1 ). (29) 
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Here we have separated the characteristic transverse dimension of the projectile 

Ri= j d 2 rr Pl (0,r) (30) 

(and a factor tt for convenience) and introduced a dimensionless function which is a 

vertex for the interaction of the projectile with a pomeron at momentum transfer y/—t± = I. 
F\(ti) behaves like In |ti | at small t\. The rest of the Green function is easily calculated by 
the stationary point method to finally give 

Xl ( s ,l,q) = -———-R 1 F 1 (t 1 )J(l,q). (31) 
yearns)' 3 / 2 

Evidently the function X2(s,l,q) is given by the same formula with the projectile form- 
factor Fx substituted by the target one F 2 . 

Combining our results for X an d Xi,2, we obtain our final expression for the absorptive 
part D corresponding to the central diffraction from Eq. (12) to be 

D = 87r- n /y 6 iV 4 (iV 2 " l)(s 2 /M 2 ) - [SlS2M > 

Vain M 2 (a In si) 3 (a In s 2 ) 3 

[RiR 2 F 1 (t 1 )F 2 (t 2 )} 2 B 2 /Vhh- (32) 

Here the number B is defined as a result of the q integration 

B = l J(d 2 q/(2ir) 2 )\l/2 + q\[V q J(l,q)] 2 . (33) 

It does not depend on I and can be represented as an integral over three momenta 

r (M(o ^4/^ f J2 ,2 ,2 l \l/2 + Q\(q + P)(Q+P')(l+P+-P-)(l+P'+-p'-) ^ 
B = l/(2vr) I) d qd pd p , , 34 

J p+p-p+P-\q + p\\q + p\ 

where 

p± = \p±l/2\, p' ± = \p' ±l/2\. (35) 
It is a well-defined integral. Calculations give 

B = 1.962 ±0.002. 

The remaining problem is to separate the PP cross-section from the absorptive part D. 

3 The pomeron-pomeron cross-section 

To separate the PP cross-section from the absorptive part for central diffractive events we 
compare the latter with the elastic amplitude A(s, t) for the scattering of the projectile and 
target in the same approach. In the single pomeron exchange approximation we have 

A(s, t) = ig A {N 2 - 1) 8 [ d 2 rd 2 r' Pl {l, r)p 2 {-l, r')G, (a, r, r'). (36) 
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The normalization is chosen to have the elastic cross-section 

da el /dt = (l/16vrs 2 )|A(s,t)| 2 . (37) 

Performing calculations similar to which lead us to Eq. (32) we find for the forward 
scattering t = 

A(s, 0) = ij^9\N- - l> 1+A RiR2 ] f^ s (38) 

and for t < 

A(s,t) = ij^9\N* - l) S ^^^ Fi(t)F2(t) _^|_ . (39) 

Evidently one cannot determine the pomeron and its couplings in a unique way from these 
expressions. One also notices that the amplitude at t = rises with s faster that at t < by 
a factor Ins. This shows that the amplitude has a discontinuous behaviour at t = 0, which 
results in the logarithmic singularity of the form-factors Fi^t) at t = 0. As a consequence 
we meet with an additional difficulty in fixing the "flux factor", since working at t < 0, we 
cannot fix the coupling by its value at t = 0. 

Therefore rather to fix the flux factor from the coupling, we choose a particular expression 
for the pomeron propagator, which simplifies the factor entering the triple pomeron vertex. 
We take that the pomeron contribution is at t = 

p( s ,o) = 2^4=r= ( 4 °) 

yams 

and at t < 

S l+A 

P ^-^J^W (41) 

Then we find for its coupling 7 to a real external particle at t = 

7(0) = 7^-5 ViV 2 - 1 R (42) 

and at t < 

l(t) = ^-j=g 2 VN^~lRF(t). (43) 

These expression define the flux factor in our case. Of course, they are not unique: one can 
always multiply the coupling by a constant, simultaneously dividing the propagator by its 
square. 

Rewriting (32) in terms of these quantities we obtain 

D = M 2 1 l{t l ) 1 l{t 2 )P 2 {siM)P 2 {s 2 M)<y%{M 2 M,t2), (44) 

where the PP total cross-section, defined by this equation, is 

^■"■"'-^^TTmP^s- (45) 
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Evidently it has the form (1) indicated in the Introduction on general grounds, with the 
constant factor c given by 

_/V 4 

c = 512^^2— j-# 2 ^ 35370. (46) 

Note that the PP cross-section can be written in terms of the triple-pomeron vertex 
l3p(t) found in [1] to be (with the flux factor given by Eqs. (40)-(43)) 

g 4 N 2 B 1 , . 

w(4) = iWTT^r (47) 

In terms of 73P and the pomeron propagator at t = (40) the PP cross-section (45) reads 

af P (M 2 MM) = ~f 3 p(hh3p(t2)P(M 2 ,0)/M 2 . (48) 

Thus the pomeron-pomeron interaction correctly factorizes into a pair of triple pomeron 
vertices joined by a pomeron which propagates between them (this fact is independent of the 



normalization of the pomeron propagators in Eqs. (40) and (41)). The characteristic I/VM2 
behaviour of the cross-section follows that of the two triple pomeron vertices. 

4 Discussion 

The obtained result corresponds to an exact prediction of the hard pomeron theory, within 
the mentioned uncertainties associated with the choice of scale in logarithmic factors. To 
relate these results to observable data one has to specify the range of transferred momenta 
to which they may apply and the value of the fixed coupling constant a s . The transferred 
momenta should not evidently be too low, to avoid confinement effects neglected in the hard 
pomeron theory. So one should take at least | *x ,2 1 > 1 (GeV/c) 2 and higher values of \t\^\ 
would be preferable. 

As to the coupling constant a s , it may naively be identified with the expected value of 
the running QCD coupling constant at the appropriate scale defined by the values of | *i,2 1 - 
However assuming that olqcd is around 0.2 at 1 GeV/c, one arrives at the well-known value of 
the intercept A ~ 0.5, which not only contradicts soft scattering data but is also about twice 
larger that the intercept observed in low -x DIS data, performed at comparable values of the 
momentum transfer [8]. For this reason we prefer to treat a s as a free parameter inherent 
in the hard pomeron theory, whose values can be deduced from the observed intercepts 
determined at appropriate momentum transfers. Low- a; DIS experiments at relatively small 
momentum transfers give the value of the intercept A ~ 0.19 [8]. From this, using (18) we 
deduce 

a s = 0.072. (49) 
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With this value of a s and at (ii^l = 1 (GeV/c) 2 we get for the triple pomeron vertex (Eq. 
(47)) 

73P (|t| = l (GeV/c) 2 ) = 0.32 mb 1/2 (50) 
and for the PP total cross-section 

0.19 

<?pp{s, 1*1.2 1 = 1 (GeV/c) 2 ) = 0.34^= mb. (51) 

yln s 

In the latter formula we assume that the appropriate scale is around 1 GeV, so that s should 
be measured in GeV 2 . At yfs = 6 TeV from (51) we get the PP cross-section 2.2 mb. 

Comparing to the existing phenomenological estimates, based on diffractive scattering 
data, we observe that the value of the triple pomeron vertex (50) is quite close to the one 
obtained from the lower energy and t data 73P = 0.364 ± 0.025 mb 1 / 2 [9]. The experimental 
value actually refers to quite small \t\ ~ 0.05 (GeV/c) 2 . Our theoretical value obtained 
for much higher |t| is only slightly smaller, which seem to indicate that the triple pomeron 
interaction freezes at \t\ of the order of 1 (GeV/c) 2 , where the confinement effects are expected 
to set in. 

As a result, our predictions of the PP cross-section are of the same order as the estimates 
done in [10] on the basis of the experimental value for the triple pomeron vertex. In this 
reference values of app for the PP cm. energy intervals 2 180 and 2 600 GeV were 
estimated to be 0.3 mb and 0.6 mb, respectively. Our Eq. (51) gives app = 0.37 mb for 
y/s = 2 GeV and 1.08 mb for ^ = 600 GeV. 

One should have in mind that the magnitude of the PP cross-section is very sensitive to 
the choice of a s and the resulting values of the intercept A. Taking a = 0.1, only slightly 
larger than (49), raises the intercept to A = 0.265, with a consequence that the PP cross- 
section at y/s = 6 TeV and \ti^\ = 1 (GeV/c) 2 rises more than 10 times to achieve values 
around 30 mb! For the range of energies considered in [10] we would then obtain PP cross- 
sections between 1.3 mb (at 2 GeV) and 9 mb (at 600 GeV). Therefore the experimental study 
of the PP cross-section presents a very efficient and sensitive test for possible signatures of 
the hard pomeron physics. 
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Figure captions 

Fig. 1. Diagram corresponding to the absorptive part for central diffractive events in 
hadronic (photonic) interactions. Ladders show pomerons emitted by the projectile and tar- 
get. 



Fig. 2. The central part of the diagram of Fig. 1. The upper and lower pomerons are joined 
to the central one via two Bartels' vertices. 
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